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Preface 


The  purpose  of  'his  thesis  was  to  study  the  effects  of  variable 
inner  wall  temperature  conditions  in  triggering  or  suppressing  complex 
multicellular  flow  fields  within  a  narrow  horizontal  annulus.  Such 
research  has  practical  applications  in  technological  systems  that  employ 
the  annular  geometry  to  cool  or  insulate  critical  components.  Some 
examples  include  nuclear  reactor  design,  materials  processing,  aircraft 
cabin  insulation  and  pressui ized-gas  underground  electric  transmission 
cables. 

In  all,  ten  inner  cylinder  temperature  distributions  were  examined, 

i 

offering  insight  into  the  nature  of  complex  annular  convective  flow 
patterns  and  the  thermal  conditions  that  cause  them.  Since  this  research 
dealt  only  with  low  Prandt 1-number  fluids,  follow-on  studies  of  variable 
Prandtl-number  effects  should  be  undertaken  to  further  explore  this 
issue . 

In  preparing  this  thesis,  I  am  deeply  grateful  for  the  assistance  of 
others  who  helped  see  me  through  this  project.  I  wish  to  thank  my 

faculty  advisor,  Capt  Daniel  B.  Fant,  whose  superb  technical  expertise  in 

/ 

this  area  filled  many  gaps  and  whose  dedicated  assistance  made  this 
effort  a  valuable  learning  experience.  I  al3o  wish  to  thank  the  computer 
science  technicians  of  Systems  Research  Laboratories  for  answering  all  my 
questions  about  the  ISC  computer  system.  Finally,  I  want  to  thank  my 
wife1MBfor  ^er  understanding  and  support  during  all  those  long  days 
and  nights  I  had  to  spend  on  the  computer. 

David  L.  Bennett 
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Abstract 


The  purpose  of  this  study  was  to  examine  natural  convection  within  a 
narrow  horizontal  annulus  subject  to  variable  inner  cylinder  temperature 
distributions.  Both  numerical  and  analytical  approaches  were  taken  in 
determining  the  effects  of  variable  temperature  conditions  imposed  on  the 
inner  cylinder  in  triggering  or  suppressing  multicellular  flow 
instabilities. 

The  two-dimensional  Navier-Stokes  equations  are  simplified  into 
boundary-layer  equations  for  the  assumptions  of  large  Rayleigh  number, 
small  annular  gap,  and  small  Prandtl  number.  These  2-D  unsteady 
boundary-layer  equations  are  discretized  using  finite-differencing 
techniques.  Numerical  solutions  to  these  governing  equations  are  then 
obtained  by  using  a  stable  second-order,  ful ly-impl ici t ,  time-accurate, 
Causs-Siedel  iterative  procedure. 

In  addition,  analytical  steady-state  solutions  to  these  simplified 
equations  are  obtained  using  perturbation  methods.  For  most  inner 
cylinder  temperature  distributions  considered  in  the  steady  flow  regime, 
these  analytical  results  yield  excellent  agreement  with  numerical 
results.  Although  both  schemes  predict  the  formation  of  thermal-like 
instabilities  induced  by  localized  hot  and  cold  spots  in  the  horizontal 
annular  regions,  the  analytical  model  failed  to  predict  the  steady-state 
multicellular  hydrodynamic  instabilities  found  numerically  for  the  sin 
temperature  distribution  at  G  =  A. 95  in  the  vertical  portions  of  the 


annulus.  The  analytical  model  also  fails  to  capture  unsteady 
■ulticellullar  flow  behavior  found  numerically  for  the  sinusoidal 
temperature  distribution  at  G  =  4.99.  Limited  unsteady  numerical  results 
for  this  test  case  indicated  oscillatory  behavior  in  the  strength  of  this 
flow  field. 
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VARIABLE  WALL  TEMPERATURE  EFFECTS  ON  MULTICELLULAR  NATURAL 
CONVECTION  IN  A  HORIZONTAL  ANNULUS 

1.  Introduction 

The  study  of  natural  convection  in  concentric  cylindrical  annuli 
receives  considerable  attention  from  researchers  due  to  its  practical 
applications  in  the  areas  of  nuclear  reactor  design,  aircraft  cabin 
insulation,  cooling  of  electronic  equipment,  thermal  storage  systems,  and 
material  processing.  The  horizontal  concentric  cylinder  geometry  is  used 
in  pressuri zed-gas  underground  electric  transmission  cables  (Kuehn  and 
Coldstein,  1975). 

The  purpose  of  this  study  was  to  examine  the  effects  of  variable 
wall  temperature  conditions  on  natural  convective  flows  in  a  narrow 
horizontal  annulus.  Of  particular  interest  is  the  effects  of  variable 
inner-cylinder  wall  temperature  distributions  in  suppressing  or 
triggering  hydrodynamic  instabilities  that  give  rise  to  the  formation  of 
convective  cell  patterns  within  the  annulus.  Understanding  the  boundary 
conditions  that  give  rise  to  these  multicellular  flow  fields  is 
important,  since  convective  cells  are  the  primary  cooling  mechanism  for 
the  annular  geometry.  Such  knowledge  may  aid  in  the  design  of  efficient 
concentric-cylinder  cooling  devices  or  related  technological  systems. 

This  thesis  involves  two  approaches  in  accomplishing  the  stated 
objective. 

First,  the  2-D  Navier-Stokes  equations  are  simplified  into 
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boundary- layer  equations  for  the  conditions  of  large  Rayleigh  nunber, 
snail  annular  gap,  and  snail  Prandtl  nunber.  These  2-D  unsteady 
boundary- layer  equations  are  discretized  using  finite-differencing 
techniques.  The  governing  equations  in  discretized  fora  are  then  solved 
using  the  AFIT  VAX  8650  conputer  systen  to  obtain  the  flow  field. 

The  second  approach  involves  deternlning  analytical  steady-state 
solutions  to  these  sinplified  equations  using  perturbation  nethods. 

These  solutions  serve  to  support  and  conpare  related  nunerical  results. 

Before  beginning  a  detailed  discussion  of  these  topics,  a  review  of 
the  related  literature  is  presented  to  exanine  recent  work  in  this  area 
of  study.  Chapters  II  -  IV  discuss  the  background  theory  and  analysis  of 
this  problea,  and  Chapters  V  and  VI  present  the  key  results  and 
conclusions. 

Literature  Review 

The  natural  convection  of  flow  between  horizontal  concentric 
cylinders  has  been  the  subject  of  aany  investigations.  In  most  of  these, 
the  focus  has  been  on  nunerical  studies  of  fluids  with  low  Rayleigh 
nuabers  or  high  Prandtl  nunbers  of  order  1  or  larger  that  show  close 
agreeaent  with  experimental  work.  The  convective  behavior  of  low 
Prandtl  number  fluids  such  as  liquid  metals  has  received  little 
attention.  These  fluids  are  often  used  as  the  working  fluid  in  several 
powergenerating  cycles  due  to  their  excellent  heat-transfer 
characteristics. 

This  thesis  considers  natural  convection  of  high  Rayleigh  nuaber/low 
Prandtl  number  fluids  in  a  narrow  horizontal  annulus  with  variable 
temperature  distributions  imposed  on  the  inner  cylinder  wall.  By 
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contrast,  most  recent  work  deals  with  isothermal  boundary  conditions, 
with  numerical  work  incorporating  variable  tenperature  effects  nearly 
nonexistent.  Analytical  approaches  in  modeling  this  flow  regime  are 
equally  rare. 

This  literature  review  covers  three  sain  topics;  the  analytical 
approaches,  the  numerical  studies,  and  the  experimental  analyses. 

Analytical  Studies.  Only  two  analytical  studies  of  natural 
convection  of  low  Prandtl -number  fluids  between  concentric  cylinders  were 
found.  Both  used  perturbation  methods  to  model  high  Rayleigh  number  flow 
fields  within  the  annulus. 

Using  the  assumption  of  small  annular  gap  width,  Fant  (1987), 
derived  a  three-term  perturbation  expansion  to  obtain  steady-state 
solutions  for  the  high  Rayleigh  number/low  Prandtl  number  flow  regime. 

He  obtained  stream  function  and  vorticity  solutions  that  agreed  well  with 
pretransitional  numerical  flow-field  results.  Since  this  thesis  adopts 
both  his  analytical  and  numerical  methods,  discussion  is  deferred  to  the 
next  three  chapters  for  a  detailed  treatment  of  this  procedure. 

Custer  and  Shaughnessy  (1977)  studied  natural  convection  of  low 
Prandtl -number  fluids  within  a  horizontal  annulus  by  solving  the 
dimensionless  thermal  energy  and  vorticity  equations  with  a  double 
perturbation  expansion  in  powers  of  Grashof  and  Prandtl  numbers.  The 
temperature  and  stream  function  expansions  were  as  follows: 
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The  authors  described  the  evolution  of  the  flow  for  several  radius 
ratios,  where  in  both  cases  the  flow  was  composed  of  a  single  eddy  in 
each  half-annulus  for  low  Crashof  number.  As  Crashof  number  was 
increased,  they  observed  that  the  eddy  in  high  Prandtl -nuaber  fluids 
fell,  while  the  eddy  in  low  Prandtl  number  fluids  rose.  Custer  and 
Shaughnessy  indicated  that  only  further  experimental  and  numerical 
studies  of  the  unsteady  flow  regime  could  resolve  multicellular  flow 
behavior. 

Numerical  Studies.  As  previously  stated,  the  majority  of  the 
numerical  studies  of  natural  convection  in  horizontal  annuli  have  dealt 
with  low  Rayleigh-nuaber/high  Prandl-nuaber  flows  under  isothermal 
boundary  conditions.  Special  emphasis  shall  be  placed  on  the  formation 
and  nature  of  multicellular  flow  fields  and  the  numerical  schemes  which 
were  used  to  model  them. 

Charrier-MoJtabi  e£.  a]_.  (1978)  studied  natural  convection  flows 
between  two  horizontal  concentric  annuli  with  constant  surface 
temperatures.  They  solved  the  problem  using  an  implicit  alternating 

direction  scheme  with  a  vorticity  and  stream  function  formulation. 

2  a 

Rayleigh  numbers  varying  between  10  and  5  x  10  ,  radius  ratios  between 
1.2  and  5  and  Prandtl  numbers  of  0.7  and  0.0Z  were  examined.  For  R  =  1.2 
and  Rayleigh  numbers  of  6000  and  7000,  no  multicellular  flow  fields 
resulted  using  air  as  the  fluid  (Pr  >  0.7).  However,  for  Pr  ■  0.02  and 
R  *  2,  the  flow  regime  became  multicellular  for  large  Rayleigh  numbers 
and  began  to  oscillate,  as  the  flow  became  unsteady. 
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Rao  at.  si-  0  985)  studied  flow  patterns  of  natural  convection  in 
horizontal  cylindrical  annuli.  They  used  an  unsteady  formulation  of  the 
2-D  energy  and  vorticity-stream  function  equations  which  were  solved  by 
using  an  ADI  scheme  with  central-differencing.  A  radius  ratio  of  1.175 
and  a  Rayleigh  number  of  4,000  resulted  in  the  formation  of  two  small 
counter-rotating  cells  just  above  the  main  kidney-shaped  cells  in  the 
upper  annulus.  Their  results  contrast  with  those  obtained  by 
Charrier-Hojtabi  et  al_.  using  the  ADI  method,  who  failed  to  achieve 
multicellular  flow  patterns  in  their  study  of  narrow-gap  geometries. 

Date  (1985)  studied  natural  convection  in  a  horizontal  annular  gap 
in  which  the  isothermal  inner  cylinder  was  hotter  than  the  cool  outer 
cylinder.  His  numerical  prediction  scheme  involved  solving  the  governing 
equations  for  the  primitive  variables  of  pressure,  velocities  and 
temperatures.  This  is  in  contrast  to  most  other  numerical  schemes  which 
use  the  stream  function-vorticity  equations  to  generate  solutions.  His 
numerical  predictions  for  natural  convection  heat  transfer  for  a  gap 
width/internal  diameter  (L/Dj)  ratio  of  0.15  for  air  compare  poorly  with 
computed  heat  transfer  rates  of  Boyd  (1981)  and  Crigull  and  Hauf  (1966). 

Farouk  and  Cuceri  (1981)  examined  both  laminar  and  turbulent  natural 
convection  between  horizontal  concentric  cylinders  held  at  constant 
temperatures.  Solutions  for  the  laminar  flow  case  were  obtained  up  to  a 
Rayleigh  number  (based  on  gap  width)  of  105.  Turbulent  flow  results  were 
obtained  for  Rayleigh  numbers  of  10®  -  107.  For  the  laminar  case,  the 
2-D  elliptic  partial  differential  equations  of  vorticity,  stream 
function,  and  temperature  with  specified  boundary  conditions  were  solved. 
An  outside-to-inside  diameter  ratio  of  2.6  was  used  with  air  in  the 
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numerical  formulation.  For  these  conditions,  the  typical  crescent-shaped 
cellular  pattern  emerged  which  has  been  found  in  previous  studies. 

Similar  fluid  circulation  patterns  were  obtained  in  the  turbulent 
studies.  However,  the  location  of  maximum  stream  function  was  found  to 
move  upward. 

Kumar  (1988)  investigated  natural  convection  in  horizontal  annuli 
for  a  wide  range  of  Rayleigh  numbers  extending  from  conduction  to  the 
convection-dominated  steady  flow  regime  for  diameter  ratios  of  1.2  -  10. 
Here,  the  inner  cylinder  is  heated  by  a  constant  heat  flux  and  the  outer 
cylinder  is  isothermally  cooled.  As  found  by  other  researchers,  a 
crescent-shaped  eddy  dominated  for  the  small  diameter  ratioo,  while  a 
kidney-shaped  flow  pattern  emerged  for  a  diameter  ratio  of  5.  At  high 
Rayleigh  numbers,  as  the  diameter  ratio  is  increased,  more  stagnant  fluid 
exists  at  the  bottom  of  the  annulus  coincident  with  the  prescence  of  the 
kidney-shaped  flow  pattern. 

Fusegi  and  Farouk  (1986)  conducted  a  three-dimensional  analysis  for 
natural  convective  flows  of  air  in  an  annulus  having  differentially 
heated  concentric  isothermal  horizontal  cylinders  using  the 
vorticity-velocity  formulation.  The  Rayleigh  numbers  used  were  103  and 
10  .  The  ratio  of  the  gap  width  to  the  inner  cylinder  radius  studied  was 
1.6.  The  authors  pointed  out  that  3-D  numerical  solutions  were  better 
suited  for  modeling  convective  flow  fields  since  2-D  numerical  studies 
could  not  capture  all  the  important  flow  characteristics  of  a  confined 
fluid.  Flow  field  results  were  graphically  compared  for  both  cases. 
Radial  velocity  plots  at  £  *  180*  indicated  that  the  3-D  flow  motion  is 
much  more  vigorous  than  the  2-D  case  due  to  the  formation  of  a 
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three-diaensional  pluae  in  this  annular  region.  This  was  but  one 
phenoaenon  Fusegi  and  Farouk  contended  that  the  2-D  solutions  overlooked. 

Xuaar  and  Keyhani  (1988)  obtained  detailed  nuaerical  results  of 
natural  convection  in  a  horizontal  annulus  with  a  constant  heat  flux 
enforced  on  the  inner  cylinder  and  an  isotheraal  outer  wall.  Prandtl 
nuabers  studied  were  0.7,  5.0  and  100  for  radius  ratios  of  1,8  <  R  <  15. 
The  Navier-Stokes  equations  were  recast  as  vorticity-streaa  function 
equations.  The  vorticity  transport  equation  and  the  energy  equation  were 
solved  using  the  false  transient  Alternating  Direction  Iaplicit  (ADI) 
aethod,  and  the  streaa  function  equation  was  solved  by  the  successive 
over-relaxation  (SOS)  aethod.  Prandtl  nuaber  effects  on  flow  field 
characteristics  in  horizontal  annuli  were  found  to  be  insignificant  for 
Pr  >  0.7  for  the  range  of  radius  ratios  studied.  The  results  also 
indicated  that  the  laportant  flow  paraaeter  for  heat  transfer  due  to 
convection  only  is  the  conventional  Rayleigh  nuaber. 

Bxperiaental  Studies.  Although  aany  articles  on  experiaental 
studies  of  natural  convection  within  enclosures  exist  in  the  literature, 
only  those  that  deal  with  the  concentric  cylinder  geoaetry  are  discussed 
below.  As  with  nuaerical  studies  of  this  problea,  experiaental  work  has 
dealt  aainly  with  the  isotheraal  condition  iaposed  on  the  annulus,  while 
variable  inner  cylinder  wall  teaperature  conditions  are  scarce  or 
nonexistent. 

Kuehn  and  Goldstein  (1976)  perforaed  an  experiaental  and  nuaerical 
study  of  natural  convection  within  horizontal  concentric  cylinders.  The 
experiaental  results  were  obtained  using  a  Hach-Zehnder  interferoaeter  to 
visualize  the  teaperature  field  in  the  enclosure  which  can  be  analyzed  to 
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determine  the  local  heat  flux.  The  flow  field  is  assumed  steady  with  no 
axial  variation  in  properties,  except  for  end  effects.  The  cylinders  had 
an  annular  gap  width/inner  cylinder  diameter  of  0.8.  Two  sets  of 
experimental  runs  were  made;  one  with  air,  the  other  using  water. 

Quantities  obtained  experimentally  included  temperature 
distributions  and  local  and  averaged  heat-transfer  coefficients.  The 
numerical  solutions  confirmed  these  experimental  results,  and  provided 
the  related  velocity  distributions  and  extended  the  results  to  lower 
Rayleigh  numbers.  The  temperature  distributions  for  both  air  and  water 
were  nearly  identical  at  similar  Rayleigh  numbers.  The  flow  was  steady 
for  all  Rayleigh  numbers  investigated  (Ra  <  105). 

Kuehn  and  Goldstein  (1976)  later  examined  the  effects  of 
eccentricity  and  Rayleigh  number  on  natural  convection  heat  transfer 
through  a  horizontal  cylindrical  annulus.  They  found  that  eccentricity 
of  the  inner  cylinder  caused  large  changes  in  the  local  heat  transfer  on 
both  cylinders,  but  the  overall  heat  transfer  coefficient  for  an 
eccentric  geometry  is  within  10  percent  of  that  for  the  concentric  case 
at  the  same  Rayleigh  number.  At  large  Rayleigh  numbers,  the  flow  was 
observed  to  become  unsteady  first  in  the  plume  above  the  inner  cylinder. 
This  unsteadiness  was  seen  to  increase  as  the  Rayleigh  number  was 
increased.  This  turbulence  rapidly  decayed  as  the  flow  moved  downward 
along  the  outer  cylinder  such  that  the  bottom  half  of  the  annulus 
remained  virtually  steady.  Hence,  laminar  and  turbulent  flows  were  seen 
to  exist  in  the  annulus  at  the  same  time. 

Sun  and  Zhang  (1986)  made  an  experimental  study  of  natural 
convection  heat  transfer  in  concentric  and  eccentric  horizontal 


cylindrical  enclosures.  They  employed  a  Mach-Zehnder  interferometer  to 
gauge  the  temperature  distribution  within  the  annulus.  The  objective  of 
this  study  was  to  determine  the  effect  of  the  diameter  ratio  K  and  the 
eccentricity  c  on  the  heat  transfer  in  the  horizontal  cylindrical 
annulus.  Experiments  were  performed  using  air  at  atmospheric  pressure 
with  radius  ratios  K  varying  from  1.77  to  2.68  and  eccentricities  from 
-0.50  to  0.35. 

The  local  and  overall  heat  transfer  coefficients  were  obtained  by 

analyzing  the  interferograms  over  a  range  of  Crashof  numbers  from 
2  4 

4.6  x  10  to  3.8  x  10  .  The  shape  of  the  interferometer  fringes 
indicated  that  the  thermal  boundary  layer  for  air  formed  around  the  inner 
cylinder  and  partly  at  the  upper  surface  of  the  outer  cylinders.  For  the 
case  of  £  -  -0.5,  temperature  inversion  and  thermal  boundary-layer 
separation  effects  were  strongest  in  the  upper  portion  of  the  annuli,  and 
to  a  leaser  degree  for  *r  =  o  and  £  =  0.353.  The  stagnation  zone  at  the 
annulus  bottom  decreased  when  e  increased  from  negative  to  positive 
values.  The  thermal  plume  formed  over  the  top  of  the  inner  cylinder 
became  more  prominent  with  increasing  eccentricity,  since  there  was  more 
space  and  less  flow  resistance  in  the  upper  portion  of  the  annulus. 

Rao  et.  §1.  (1985)  experimentally  obtained  results  for  steady  3-D 
natural  convection  in  a  horizontal  cylindrical  annulus.  The  objective 
was  to  confirm  the  complex  structure  of  the  3-D  numerical  results 
obtained  by  the  authors  of  this  paper.  Flow  patterns  were  examined  for 
the  large  Rayleigh  number/large  Prandtl  number  regime.  An  increase  in 
the  number  of  cells  observed  with  increasing  Ra  confirmed  numerical 
studies  of  similar  flows. 
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II .  Mathematical  Analysis 


In  this  chapter,  a  vorticlty-stream  function  formulation  of  the 
Navier-Stokes  equations  is  described  in  dimensionless  form  (based  on 
Fant,  1987).  By  using  this  approach,  the  pressure  terms  are  completely 
eliminated  while  automatically  satisfying  the  conservation  of  mass 
principle.  In  addition,  the  resulting  stream  function  and  vorticity 
contours  are  useful  for  visualizing  and  analyzing  the  flow  field.  The 
dimensionless  equations  in  final  form  are  given  in  Eqs  (2.26  to  2.28). 


The  Physical  Model 

a.  The  flow  is  unsteady  and  two-dimensional  (see  Figure  2.1). 

b.  The  fluid  is  initially  at  rest. 

c.  The  cylinders  are  assumed  horizontal  with  an  Isothermal  outer 
cylinder  and  a  variable  temperature  distribution  imposed  on  the 
inner  cylinder. 

d.  Laminar  fluid  motion  is  induced  by  buoyancy  effects.  The  fluid 
is  Newtonian. 

e.  All  material  properties  are  assumed  constant.  Density 
variations  are  allowed  to  occur  via  the  Boussinesq 
approximation. 


The  Dimensional  Formulation 

Governing  Equations  Using  the  Stream  Function-Vortlcitv  Approach. 
The  vorticlty-transport  equation  is  derived  by  taking  the  curl  of  the 
Navier-Stokes  equations.  This  equation  can  be  written  as: 
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(2.1) 


-  -  -  -  2--- 
—  ♦  (v  •  7)w  =  (w  •  7)v  ♦  vV  w  4  7  x  F 

at 

The  bar  over  Che  variables  signifies  dimensional  quantities,  w  and  v 

represent  the  vorticity  and  velocity  vectors,  respectively.  The  left 

side  of  this  equation  desribes  the  total  rate  change  of  particle 

vorticity.  The  ter*  (w  •  7)v  is  identically  zero  for  2-D  flows  due  to 

the  vorticity  vector  always  being  perpendicular  to  the  plane  of  flow. 

2- 

The  tere  t>7  w  represents  the  net  rate  of  vorticity  diffusion  due  to 
viscous  effects.  The  last  tere,  7  x  Fg,  is  the  rate  of  internal 
vorticity  generation  due  to  body  forces,  resulting  from  large  densi  ty 
perturbations  in  natural  convective  flows.  The  fact  that  pressure  does 
not  appear  explicitly  in  Eq  (2.1)  enables  the  deternination  of  the 
vorticity  and  streae  function  (velocity)  fields  without  prior  knowledge 
of  the  pressure  distribution.  Using  the  Boussinesq  approximation  for  the 

net  body  force: 


-  Tq)T  =  -gtf(T  -  Tq)cos  v  e r  4  g/?(T  -  TQ)sin  W  e^  (2.2) 

where  ft  is  the  coefficient  of  thermal  expansion.  The  final  form  of  the 
vorticity-transport  equation  is: 

5W  1  d(f,w)  ,  &T  cos  V  0T  ^  j- 

—  *  -  — 7-7-  =  gft  I  sin  V  —  +  — ^ - 7  I  +  vV  w  (2.3) 

0t  r  d(r,y)  r  * 

and  the  thermal -energy  and  stream  function  equations  are: 


and 


0T  1  0(f,T) 

—  + - -  aVT 

Ot  r  0(r,y) 

w  =  T2!  . 


(2.4) 

(2.5) 
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Ia  the  above  equations,  r  and  represent  the  radial  and  angular 
coordinates.  T  is  the  temperature  and  t  denotes  time,  v  and  a  denote 
the  momentum  and  thermal  diffusivities,  respectively.  The  radial 
velocity,  u,  and  the  angular  velocity,  v,  have  been  written  in  terms  of 
the  stream  function,  f,  defined  as: 


-1  9f  9f 

u  - - and  v  =  — 


(2.6) 


r  dys  dr 

The  number  of  dependent  variables  has  now  been  reduced  by  one.  The 
convective  terms  are  then  given  in  terms  of  the  Jacobian,  9(P,Q)/9(x,y) , 
where 


*<P,Q)  9?  9Q  9?  9Q 

-  - - .  (2.7) 


*<x,y)  9x  9y  9y  9* 

The  three  coupled  governing  equations  (2.3,  2.4  and  2.5)  describe 
the  vorticity,  temperature  and  stream  function  for  the  horizontal 
annulus. 

Boundary  Conditions.  For  temperature  boundary  conditions,  the  outer 
cylinder  is  assumed  to  be  isothermal,  while  the  inner  cylinder  has  a 
variable  temperature  distribution  that  is  a  function  of  the  angular 
position. 


T  =  Tj(yo  at  r  *  a 

*nd  T  »  Tfl  at  r  =  b  .  (2.8) 

Both  u  and  v  are  identically  zero  at  both  cylinder  walls  due  to  the 
no-slip  condition,  yielding  the  boundary  conditions: 
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1  =  0,  w 


— =-  at  r  =  a  and  r  =  b. 
Or1 


(2.9) 


To  ensure  a  computational  continuous  condition  in  the  complete  annulus 
(0  -  In)  the  following  boundary  conditions  (in  yt)  are  enforced  at  zero 
and  In  radians: 


and 


for  a  £  r  <  b, 


(2.10) 


*  B  *2n 
9<p  0<p 

—  =  T  I 

9y>  0  1 2n 

where  <p  represents  the  dependent  variable  for  temperature,  vorticity,  or 
stream  function. 

Initial  Conditions.  Since  there  is  no  fluid  motion  initially, 


f  =  w  =  0  and  T  =  T0  (2.11) 

throughout  the  annulus,  except  at  the  walls  where: 

T  *  Tj(v)  and  T  =  Tq  (2.12) 

is  enforced  at  r  =  a  and  b,  respectively. 


The  Oimensionless  Formulation 

Casting  equations  in  dimensionless  form  has  several  advantages. 
First,  it  transforms  the  mathematical  or  numerical  results  into  a  simpler 
form,  allowing  for  better  graphical  interpretation.  Secondly, 
measurement  scales  are  no  longer  an  intrinsic  part  of  the  physical 
quantities,  therefore,  any  laws  governing  physical  variables  are  valid 
for  all  different  measurement  systems  (Fant,  1987).  Finally,  when  a 
problem  is  nond i mens ionali zed,  fewer  variables  are  used  and  the  proper 
dimensionless  groups  characterizing  a  particular  flow  model  are  usually 
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brought  forth. 


Coordinate  Transformation.  The  radial  coordinate  is 
nond i mens ionali zed  so  that  the  outer  boundary  at  r  *  b  is  transformed 
into  the  unit  circle,  r  =  1.  The  inner  boundary,  r  =  a,  is  transformed 
into  the  pole,  r  =  0.  The  following  coordinate  transformation  results: 

r  -  a 

r  ■  -  .  ?  =  ?  .  (2.13) 

b  -  a 

The  other  independent  variable,  time,  scales  as 


tL> 

t  =  —  .  (2.14) 

& 

Governing  Equations.  The  remaining  dependent  variables  can  now  be 
scaled  to  produce  the  following: 


T  = 


V*>  -  T0 


(temperature) 


(vorticity) 


(2.15) 


(2.16) 


f 

f  =  —  (stream  function)  (2.17) 

v 

Eqs  (2.13  to  2.17)  are  now  substituted  into  Eqs  (2.3  to  2.5)  resulting  in 
the  following  nondimensional  governing  equations: 


(2.18) 
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&I 

,  1  ^-1  <»(f,T)  1 

f  *2T  1 

—  + 

f  r  +  ] 

(2.19) 

at 

L  G  J  a<r,r)  Pr 

l  2  J 

where 


(2.20) 


Also, 

v 

Pr  =  -  (Prandtl  nuaber)  (2.21) 

a 


va 


(Rayleigh  nuaber 
based  on  the  inner 
radius,  a) 


(2.22) 


b  -  a 

and  G  =  -  (gap  nuaber)  .  (2.23) 

a 

The  three  diaensionless  paraaeters  that  arise  froa  this  analysis  are  G, 

Pr  and  Ra.  These  key  variables  are  used  to  siaulate  various  flow 
conditions  and  geoaetries.  The  aagnitude  of  the  Rayleigh  nuaber  reflects 
the  ratio  of  destabilizing  buoyant  forces  to  the  aore  stabilizing  viscous 
forces.  It  also  predicts  the  onset  of  thernal  and/or  hydrodynaaic 
instability.  The  Prandtl  nuaber  characterizes  the  fluid  and  the  gap 
nuaber  defines  the  annular  geoaetry. 

In  deriving  Eqs  (2.18  to  2.20),  the  variables  w,  f  and  t  are 
non-diaenslonalized  with  v  instead  of  a.  This  was  done  in  order  to 
reflect  a  theraal-diffusion  doainated  energy  equation  for  saall  Prandtl 


1 


I 
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nuabers,  which  was  the  focus  of  this  study.  For  further  details  on  this 
procedure  see  Fant  (1987). 

In  8qs  (2.18  to  2.20)  the  Jacobian  terns  can  be  expressed  In  terns 
of  the  velocity  variables,  u  and  v,  where 

1  f  1  i"1  -  f  8 

u  =  —  I  r  +  -  I  —  =  u  I  —  I  (radial  velocity)  (2.24) 

C  GJ  9? 


and 


1  f  r  a 

v  = - =  v 

G  r  "  v 


(:) 


(tangential  velocity)  (2.25) 


With  these  substitutions,  the  diaensionless  governing  equations  becoae: 
Energy: 

91  91  ,  1  ^-1  91 

G  '  9? 


C  —  +  uG 

at  Or 


91  rly 

—  +  vG  r  +  - 
jtw  l  r.  J 


+ 1 

[r.l 

.-I  «T 

1  —  + 

)'2^} 

(2.26) 

l  «r 

^  G  ■ 

*  9c 

G  ■ 

’  tv1  J 

Vorticity: 


{  Ou  9»  f  U-l  ] 
G  —  +uC  —  ♦vC|r  +  -|  —  > 

9t  9c  G  J  9f>  ) 


r  a2w 

r  1  ■ 

-i  aw  1 

»  Pr  I  — r-  +  uG  —  +  vG 

[  r  +  - 

—  \ 

[fir1  9c 

L  G  ■ 

>  9?  J 

+  G(Ra) 


r  *t 

r  1 

-1  91  1 

4  sin*'  —  4  cosy* 

r  ♦  - 

- 

(2.27) 

*■  G 

>  9p  J 
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Stress  Function: 


— 2  + 

(r.i‘ 

fr.i- 

-2  ah 

— J  =  o »  • 

(2.28) 

Or1  1 

G 

*  Or 

G  ■ 

>  Ov2 

The  fori  of  Eqs  (2.26  -  2.28)  facilitates  the  asymptotic  analysis  that 
follows  in  Chapter  III,  where  the  limits  of  Ra  -*  oo,  G  -»  0  and  Pr  -»  0  are 
separately  exaeined. 

Boundary  Conditions.  The  nondimensional  boundary  conditions  for 
this  problem  are  obtained  by  substituting  Eqs  (2.13  -  2.17)  into  Eqs 
(2.8  -  2.10).  They  are  at  r  =  0  (inner  cylinder): 

T  = 

w  = 

f  * 

and  at  r  *  1  (outer  cylinder): 

T  = 


V*) 


1  02f 
G2  Or1 


(2.29) 


f  =  0  .  (2.30) 

The  computational  continuous  condition  at  zero  and  In  radians  still 
applies,  except  that  dimensional  quantities  are  now  nondimensional i zed. 

Initial  Conditions.  The  dimensionless  initial  conditions  are 
virtually  identical  to  Eqs  (2.11  and  2.12)  for  the  entire  annulus: 

T  -  f  «  h  -  0f  (2.31) 

except  at  the  inner  and  outer  cylinder  walls,  where: 
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and 


T  *  Tjly)  at  r  =  0 

T  =  0  at  r  =  1  .  (2.32) 

These  initial  conditions  assume  a  motionless  fluid  (Ra  *  0),  where  the 
problem  becomes  the  steady-state  conduction  case  for  a  cylindrical 
annulus.  However,  when  numerically  searching  for  transitional  or 
unsteady  flow  behavior,  initial  conditions  other  than  zero  are  used. 
These  procedures  are  discussed  in  Chapter  IV. 
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Ill .  Analytical  Analysis 

In  this  section,  simplified  boundary-layer  equations  are  derived 
from  the  2-D  Navier-Stokes  equations  for  certain  limiting  conditions  (see 
Fant,  1987).  Analytical  steady-state  solutions  to  these  equations  are 
then  determined. 


The  Boundary  Layer  Equations 

Using  the  governing  Navier-Stokes  equations  as  shown  in  Chapter  II 
(Eqs  2.26  to  2.28),  Fant  (1987)  obtained  the  following  expansions  of  the 
dependent  variables  for  the  conditions  of  infinite  Rayleigh  number 
(Ra  -»  oo)  and  small  annular  gap  width  (C  -»  0): 


u  =  Ra1/4u  +  0(Ra"l/4) 


v  =  Ra1/2v  +0(1) 


w  =  Ra3/4w  +  0(Ra1/4) 


f  =  Ral/4f  +  0(Ra"1/4) 


T  *  T  +  0(Ra'1/2)  .  (3.1) 

Using  these  expansions,  the  Navier-Stokes  equations  reduce  to  the 
following  Cartesian-like,  boundary- layer  equations: 

Energy: 

„  91  „  r  Of  91  Of  OT  }  1  02T 

C  —  +  G  < - + - \  - - =-  (3.2) 

Ot  ^  Oy  Or  Or  Oy  J  Pr  Or* 
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Vortlcity: 


dw  f  df  dw  df  dw  1  d2w  dT 

G2  •  Pr  -^  +  G  •  Pr  < - + - [  =  Pr— y  +  G  sinw  —  (3.3) 

dt  dr  dr  dy/  J  dr  dr 

Stream  Function: 


— ?  * 

C2w 

(3.4) 

dr2 

with  related  boundary  conditions, 

T(0,y)  *  TjW, 

T(1,V0  = 

T0(V) 

(3.5) 

f(0,y)  =  f(l,vO 

=  0 

(3.6) 

2~ 
l  d*f 

w(0,*/)  *  w(l,y) 

~  ~7  ~  9 

(3.7) 

C2  dr2 

r  =  0,1 

Note  that  the  annular  curvature  effects  appear  through  sin?'  in  the 
buoyancy  ter*  of  the  vorticity  equation  above.  Variable  wall  temperature 
Influence  comes  in  through  Tj  and  TQ  in  Eq  (3.5).  These  governing 
equations  are  valid  for  finite  Prandtl  number.  Importantly,  the  Rayleigh 
number  and  gap  number  dependency  has  vanished  and  is  replaced  by  a  single 
scaled  gap  parameter 

G  =  Ra1/4  C  .  (3.8) 

To  maintain  the  unsteady  term,  time  was  scaled  as 

t  =  t  Ra1/2  .  (3.9) 

The  scaling  of  these  equations  is  based  on  the  scaling  typical  of  laminar 


natural  convective  boundary- layer  flows,  which  cone  forth  under  the 
conditions  of  high  Rayleigh  number  and  snail-gap  width,  nanely, 

6  ~  Ra"1/4  .  (3.10) 

In  this  study,  the  finite  Prandtl  nuaber  boundary-layer  equations 
are  siaplified  for  the  Halting  condition  of  ?r  -*  0,  characteristic  of 
liquid-aetal  fluids.  As  discussed  in  the  literature  review,  other 
studies  of  natural  convection  in  narrow  annuli  have  aostly  concentrated 
on  air  as  the  working  fluid.  The  study  of  snail  Prandtl  nuaber  effects 
on  flow  instability  for  this  particular  geoaetry  are  rare,  especially 
with  regard  to  the  variable  wall  teaperature  influence.  The  derivation 
that  follows  addresses  this  issue. 

First,  let  vorticity  and  streaa  function  scale  as: 

w  •  X^  +  -  (3.11) 

f  =  X2F  + .  (3.12) 

The  streaa  function  equation  (3.4)  requires 

X2  ~  ^2>'l  ‘  <3.13) 

With  Eqs  (3.11)  and  (3.12),  the  vorticity  equation  (3.3)  becoaes 

0W  „  f  <9F  <JW  <JF  <?W  ) 

GZPrX  —  +  Pr  GX  X  l - + - f 

Ot  Z  1  &r  *  iV  J 

92U  „  91 

=  Pr  X.  — =•  +  G  sinv  —  (3.14) 

Ox1  9r 
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or, 


C2PrX 


1  ~ 


Pr  C  \2  \l  ~  Pr  ~  G 


(3.15) 


in  order  to  retain  the  physics  of  the  problea.  Balancing  the  viscous 
tern  with  the  buoyancy  tern  gives 


x  ~  _ 

1  — 
Pr 


then  Eq  (3.13)  yields, 


(3.16) 


X2  * 
Pr 


(3.17) 


Eqs  (3.16)  and  (3.17)  are  then  substituted  into  Eq  (3.15)  to  get 

G  ~  0(Pr1/4)  (3.18a) 

~  ~2 

and  t  ~  G  ,  which  yields  froa  Eq  (3.18a) 

t  ~  0(Prl/Z)  .  (3.18b) 

With  these  relations,  new  scaled  variables  can  be  defined: 


G  =  G  Pr 


1/4 


(3.19a) 


and 


Therefore, 


t  =  t  Pr 


1/2 


(3.19b) 


~  Pr 


-3/4 


(3.20) 
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(3.21) 


\2  ~  Pr 


•1/4 


which  yield  the  following  expansions: 


w  ~  Pr-3/4  W  ♦ 

0(Pr1/4> 

(3.22) 

f  ~  Pr‘1/A  F  ♦ 

0<Pr3/4> 

(3.23) 

T  ~  T  +  0(Pr) 

(3.24) 

Scaled  gap  nuaber  and  tiae  can  now  be  expressed  in  teras  of  the  key 
dimensionless  variables: 


C 


t 


(3.25) 


(3.26) 


These  new  scaled  variables  show  a  unique  dependency  on  Ra,  Pr,  C  and  t. 
Substitution  of  Eqs  (3.22  -  3.26)  into  Eqs  (3.2  -  3.4)  yields  the 
following  siaplified  set  of  governing  equations: 

Energy: 


(3.27) 


Vorticity: 

,2  „  f  <?F  OU  O?  0W  I  92 W  „  91 

C  —  +  G< - + - [  =  — =■  +  G  siny  —  (3.28) 

dl  9y  9t  9t  9y/  J  &vL  9t 

Stream  Function: 


— 7  =  G2  W  .  (3.29) 

9c 
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The  energy  equation  is  completely  uncoupled  from  the  vorticlty  equation. 
The  boundary  conditions  for  the  energy  equation  are  as  follows, 


T(0,*')  =  Tjiy) 


(3.30) 


T<  I,*')  =  Tq(v/) 


(3.31) 


where  T^(y/)  and  T^yO  denote  inner  and  outer  wall  temperatures, 
respectively.  The  isothermal  outer  wall  temperature  for  this  model  is 
taken  to  be  zero,  while  the  inner  wall  remains  variable,  or 


T(0,vO  =  Tj  ( V') 


(3.32) 


T(l,yO  =  0 


Integrating  equation  (3.27)  twice  yields, 


T  *  C^r  +  C2 


(3.33) 


(3.34) 


and  applying  the  boundary  conditions  of  Eqs  (3.32)  and  (3.33)  results  in 
the  following  expansion  for  temperature  across  the  annulus: 


and 


T(r,y)  =  (Tq  -  Tj)  r  +  Tj 


91 


(3.35) 


(3.36) 


Using  Eq  (3.36),  the  governing  equations  reduce  to  their  final  form: 
Vorticlty: 

*2,f  . 

(3.37) 


G 


<5W  „  <?F  «9W  &F  «JW  1  <?  W 

+  G  - - + - r  =  ~ 

dt  9v  9t  J  9t 


+  G  sin  y>  (T  -  T. ) 

o  1 


Stream  Function: 
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(3.38) 


a 

G  W 


with  boundary  conditions 

F(  0  ,y/)  =  F(  1 ,  V*)  =  0  (3.39) 

ff¥  m 

—  (0  ,v)  =  —  (l,v/)  =  o  .  (3.40) 

9r  &c 

These  equations  will  be  solved  nuaerically  in  Chapter  S  using  various 
inner  wall  temperature  profiles  in  order  to  study  flow  behavior  and  the 
aulticellular  instabilities  that  aay  arise. 


The  Steady-State  Perturbation  Solution 

Rearranging  Eq  (3.38)  results  in  the  following  expression  for 
vorticity: 


W  = 


(3.41) 


Assuaing  steady-state  conditions  and  substituting  Eq  (3.41)  in  Eq  (3.37) 
yields: 


{*F  *3F  *F  a3F  ’l  tf*F  . 

- -  + - - -  l  =  — -  +  (T  -  T. )  sin  w  • 

9t  9t  )  Or  °  1 


(3.42) 


Taking  the  Unit  of  C  ■*  0,  one  can  derive  a  perturbative  solution  to  the 
governing  equations.  By  exaaining  Eq  (3.42)  one  can  deduce  the  expansion 
for  streaa  function  to  be 


*3  *7  *11 

F  =  C  tj  t  C  Fj  f  0(C") 


(3.43) 
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and  fron  Eq  (3.41),  vortlclty  can  be  expressed  as 


H  «  C  Nj  ♦  C5  W2  *  0(C9)  . 


(3.44) 


The  first  effects  of  the  convective  terns  cone  in  through  and  U^. 
Substituting  expansion  (3.43)  into  Eq  (3.42),  and  Batching  terns  of  like 
coefficients  yields  the  following  steady-state  result: 


(3.45) 


*  7 

*r4 

C  : 

<To 

-  V 

*F1 

+  !!i 

Or4 

Or3 

Or 

(3.46) 


Integrating  Eq  (3.45)  yields, 

4  3  2 

r  r  r 

F.  =  —  <Tt  -  T  )  sin  y/  4  C.  —  4  C,  —  4  C,  r  4  C,  .  (3.47) 

1  4!  1  0  1  31  2  2  3  4 

The  constants  of  integration  are  obtained  fron  the  boundary  conditions  of 
Eqs  (3.39)  and  (3.40)  and  result  in: 


F,  *  r2  (r  -  1)Z  (T.  -  T  ) 
1  l  o 


sin  y 


(3.48) 


Differentiating  Eq  (3.47)  twice  with  respect  to  r  yields 


r  (r  -  1)  4 


sin  y/ 


(3.49) 


Substituting  Eq  (3.48)  Into  Eq  (3.46)  yields 
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4  3  2 

15r  +  14rJ  -  6r  ♦  r 


} (I' 


V 


2  sin  2w 
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(3.50) 


Integrating  Eq  (3.50)  four  tines  and  applying  boundary  conditions 
(Eqs  3.39  and  3.40)  yields 


F 


2 


sin  2y> 
96 


336  180 


6 

r 

180 


360  1512  5040  J 


and  becomes 

.  sin  2y/  f  r7  r6  7r*  r*  r3  r  IT 

W  =  (T  -  T  j  -  4 - ♦ - + - +  -  \  (3.52) 

96  [  21  6  30  6  18  252  2520  J 

The  expressions  for  F2  and  W2  contain  the  tern  sin  2*/,  which  nay 
contribute  to  the  multicellular  behavior  of  the  flow  field  at  some 
critical  value  of  C.  The  one-tern  expansions,  found  to  be 

-3  2  2  Sin  ^ 

F  =  GJ  r  (r  -  IT  <Tt  -  T  )  -  (3.53) 

I  o 


{1  ")  sin  V' 

r  (r  -  1)  +  -  V  (Tt  -  T  )  - , 

6  J  1  0  2 


(3.54) 


provide  good  approximations  in  modeling  pretransitional  flow  behavior  for 
Pr  -*  0.  Since  many  terns  in  the  perturbation  expansion  nay  be  required 
to  capture  complex  flow  behavior,  one  expects  the  two-tern  equations  to 
be  valid  for  only  weakly-nonlinear  flow  Instability.  In  this  study,  the 
full  effects  of  the  nonlinear  convective  terns  were  handled  by  solving 
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the  partial  diferential  equations  numerically  (see  Chapter  4). 

Eqs  (3. S3)  and  (3.S4)  represent  analytical  solutions  for  the 
steady-state  pretransitional  stream  function  and  vorticity,  applicable 
throughout  the  annulus  (0  -  In).  In  this  form,  a  variety  of  variable 
wall  temperature  conditions  can  be  considered.  These  temperature 
profiles,  which  can  vary  with  the  angular  coordinate  y,  may  be  used  to 
model  various  hot  and  cold  spots  on  the  cylinder  walls.  Several 
interesting  variations  of  the  inner  cylinder  wall  temperature  are 
described  in  detail  in  Chapter  S.  Numerical  and  analytical  solutions 
relating  to  each  inner  wall  temperature  variation  are  also  presented  and 
compared  in  this  chapter. 
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There  are  two  topics  of  discussion  in  this  chapter.  The  first 
describes  the  finite-differencing  method  used  in  solving  the  coupled  set 
of  governing  partial  differential  equations.  These  equations  are  the  2-D 
unsteady  boundary- layer  equations  developed  in  Chapter  3.  The  second 
topic  deals  with  the  computational  procedure  employed  to  solve  the 
resulting  finite-difference  equations.  Here,  the  iteration  sequence, 
convergence  criterion,  and  relaxation  parameters  for  the  numerical 
procedure  are  described. 

The  basic  numerical  approach,  relating  to  the  zero-Prandtl  number 
flow  regime,  was  adopted  from  Fant  (1987).  In  this  study,  however,  the 
buoyancy  term  in  the  vorticity  equation  was  modified  to  account  for 
nonisothermal  effects. 

In  handling  the  streamwise  nonlinear  convective  terms,  a  corrected 
second-order  upwind  scheme  was  used  for  the  boundary- layer  equations. 

For  the  radial  convective  terms,  a  corrected  second-order  central 
difference  expression  was  employed.  This  approach  ensured  numerical 
stability  when  solving  highly  convective  flow  problems.  The  unsteady 
form  of  the  equations  provided  the  opportunity  to  capture  both  steady  and 
unsteady  behavior  when  solved  in  a  time-accurate  fashion. 

The  Humeri cal  Method 

The  2-D  unsteady  boundary-layer  equations  presented  in  Eqs  (3.37) 
and  (3.38)  are  discretized  using  finite-differencing  techniques.  The 
governing  equations  in  discretized  form  are  well-suited  for  solving  on  a 


high-speed  computer  system.  The  computational  domain  consists  of  a 
cellular  mesh  formed  by  the  intersection  of  variably  spaced  radial  lines 
with  circular  arcs  concentric  with  r  *  1.  The  grid  nodes  are  the  points 
where  these  radial  lines  and  concentric  arcs  intersect.  The  cellular 
mesh  generalized  for  variable  increments  is  shown  in  Figure  4.1.  Since 
the  angular  distance  between  the  nodes  was  kept  constant  in  this  model, 
the  variables  (K^,  ...)  in  Figure  4.1  were  simply  replaced  by  the 

constant  K.  In  general,  the  use  of  variable  increments  in  the 
computational  domain  permits  the  concentration  of  nodes  in  areas  of  large 
gradients,  such  as  the  boundary- layer  regions  near  the  inner  and  outer 
cylinder  walls. 

For  all  temperature  distributions  imposed  on  the  inner  cylinder,  the 
complete  annular  flow  field  (0  -  2n)  was  numerically  resolved.  Vertical 
symmetry  was  not  assumed. 

Variable  Increment  Finite-difference  Formulas.  The  derivatives  in 
the  governing  equations  are  transformed  into  variable  increment 
finite-difference  expressions  using  Taylor  series  expansions. 

Standard  expressions  for  forward,  backward  and  central  differences, 
taking  into  account  formal  truncation  errors  for  variable  increments,  are 
given  in  Fant  (1987). 

For  the  unsteady  term  in  the  equations,  a  stable  forward-difference 
molecule  is  used: 


9$ 

a" 

«t 


0 


* 


n+1 


-  4> 

▼Vi 


♦  out) 


At 


(4.1) 


where  n  is  the  time  level  at  which  the  dependent  variable  is  evaluated. 
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gure  4.1.  Coaputational  mesh  with  variable  increments 


Finite-Difference  Equations  for  the  Dependent  Variables.  All 
spatial  derivatives  in  the  governing  equations  are  second-order  centrally 
(or  upwind)  differenced.  This  Includes  the  convective  teres  in  the 
vorticity  equation  (3.37)  which  are  represented  by  a  first-order  upwind 
expression  together  with  a  correction  term  for  second-order  accuracy. 

The  convective  teres  are  split  in  this  Banner  in  order  to  enhance 
nueerical  stability  when  resolving  secondary  flow  behavior.  This 
procedure  is  explained  in  greater  detail  in  Fant  (1987).  For  the 
zero-Prandtl  boundary-layer  equations,  only  the  vorticity  and  streae 
function  equations  (Eqs  3.37  and  3.38)  had  to  be  solved  nueerically  in  a 
coupled  Banner.  The  following  one-equation  fornat  is  used  to  represent 
these  two  coupled  equations: 

f  9*  I  *4>  *4> 

A2  <  -  AG  -r  ♦  — j  V  +  2X  —  +2*J  —  +  S  =  0  (4.2) 

[  91  9r  J  9x 

with  the  nonlinear  convective  teres  expressed  as: 


9$ 

2X  —  =  (X  -  | X | ) 

9x 


-  -  <1  *  HI 


hb  ♦  hf 


«•  (X  +  |X|) 


(— 

l  kf 


-  <PQ  ^/H  +  -  (1  +  1/H)  *0 


hb  +  hf 


) 


(4.3) 


and 
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where  4>  represents  the  dependent  variables  T  or  W,  and 


Also, 


2\  =  A1 


*F 

&141 


and 


2p  =  -  A1 


W 

9t 


where 


|  C,  <p  .  W 

1  0,  4>  =  F 


1,  <P  =  w 
1,  ^  =  F 


with 


f  1,  *  -  W 

l  0,  *  -  F  . 


(4. A) 


(4.5a) 


(4.5b) 
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For  the  radial  convective  teras  (Eq  4.3),  the  first  teras  within  the 
brackets  are  the  first-order  upwind-difference  coaponents  and  the  second 
teras  are  the  added  corrections  which  bring  the  differencing  up  to 
second-order  accuracy.  The  X.  teras  within  the  parenthesis  ensure  that 
stable  differencing  * into-the-wind'  is  enforced,  since  the 
boundary- layer  equations  are  only  first  order  in  a  corrected 
second-order  upwind-differenced  expression  was  used  to  represent  the 
streaawise  convective  terns  (Eq  4.4).  Here,  the  second  tern  within  the 
brackets  represents  the  added  correction  to  achieve  second-order 
accuracy. 

Osing  the  above  results,  the  boundary- layer  equations  are  written  in 
the  following  finite-difference  fora: 

<A2)  A  G2 

^  n+1  u.  n  .  „  ^  n4l  .  „  ^  n4l  .  „  ^  n4l 

C0*0  *  -  +  C x*i  4  C  fa  +  C 3*3 

T 


♦  c>  *  e"*1  «  sn+I 

4  4 


(4.6a) 


where 


( A2)  A  G  2(A2)  f  1  1  ^  r  1  1  ^ 

- + - +  r\ - +  riM  I  —  +  —  I 

t  hbhf  ^  hf  hfe  J  v  hb  hf  J 


4  r  |p|  £  -  J  (4.6b) 
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In  addition, 


,ntl 


f  - 

C  sin  y>  {  T 


OUTER  (VV 


-  c2  w0n+1, 


tihhbr(vo)  *  s  w 


=  f 


(4.6i) 


where  the  buoyancy  tern  (Eq  4.61)  permits  variations  in  temperature  with 
V  on  the  inner  and  outer  cylinders. 

Boundary  Conditions.  Considering  the  complete  annulus,  the 
finite-difference  form  of  the  boundary  conditions  are: 

Energy: 


(4.7a) 


Vorticity: 


Tn+1  _  A 

NR, J  " 


2Fn+1 
un+1  -  2,j 

1  i  “  '  2 

,J  (Gh^ 
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,^1+1  HR-1,  j 
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Stream  Function: 


■  • 


(4.7b) 


(4.7c) 


(4 . 7d) 


(4. 7e) 


•  (4-7f) 

To  satisfy  continuity  at  zero  and  2n  radians,  a  computational  continuous 
condition  must  be  defined  such  that: 
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where  $  represents  the  temperature,  vorticity,  or  stream  function. 

Again,  for  further  discussion  of  this  differencing  procedure,  see  Fant 
(1987). 

The  Computational  Procedure 

The  system  of  coupled  finite-difference  equations  (as  described  in 
the  preceding  section)  are  solved  implicitly  in  time  using  a  point 
iterative  Causs-Siedel  method  with  underrelaxation.  At  a  given  time 
level,  the  dependent  variables  are  found  by  repeated  iterations  of  the 
governing  equations. 

With  the  inner  and  outer  wall  temperature  distributions  preset,  the 
dependent  variables  were  initialized  by  setting  them  to  zero.  Previously 
converged  results  at  a  particular  time  level  were  used  as  initial 
conditions  to  start  the  next  time  level  calculations.  This  procedure  was 
repeated  until  steady-state  was  achieved  or  unsteady  behavior  developed. 

Iteration  Sequence  and  Convergence  Criterion.  The  boundary- layer 
governing  equations  in  finite-difference  form  were  numerically  iterated 
in  the  following  order: 

i.  vorticity  equation 
il.  stream  function  equation. 

For  the  small-Prandtl  number  limiting  condition,  the  energy  equation  was 
integrated  analytically. 

The  above  sequence  was  repeated  until  successive  iterations  were 
within  a  prescribed  tolerance.  The  maximum  modulus  of  the  difference 


between  vorticity  and  streaa  function  for  two  successive  iteration  values 
defined  the  relative  convergence  constraint,  or 


.■4-1, n+1  ,a,n+l 
<P  ~  <P 

(^B+1»n+1) 


<  1  x  10 


(4.8) 


where  4>  again  represents  vorticity  or  streaa  function.  When  Eq  (4.8)  was 
satisfied,  the  nuaerical  solution  was  considered  converged.  Here,  a 
refers  to  the  iteration  level  and  n  denotes  the  tiae  level. 

Relaxation  Paraaeters .  Two  independent  relaxation  paraaeters, 
and  were  incorporated  into  the  finite-difference  equations  to  control 
the  rate  of  convergence  of  the  iterations.  The  relaxation  paraaeter,  Oj , 
was  used  in  conjunction  with  the  vorticity  boundary  conditions.  For 
exaaple,  on  the  inner  cylinder: 


wBtl,nfl  =  ^  {  2Fa,n+l/(“h^2  }  +  (1  _  (4.9) 

along  r  =  0,  where  denotes  the  radial  spacing  between  the  first 

adjacent  node  and  the  wall.  Typically,  .1  <  fij  <  .5  was  the  range  used 

\ 

to  help  stabilize  the  nuaerical  coaputations.  The  second  relaxation 
paraaeter,  CL,  was  used  with  the  second-order  upwind  differenced 
correction  teres  which  were  based  on  the  new  tiae  level  (En+1).  The 
Interior  point  equations  for  the  vorticity  were  coaputed  according  to: 

<A2)A  G2 

„af  1  ,nf  1  ,,n  .  ,  n  ,,a,n+l  .  „  „a,n+l  .  „  „a+l,n+l  .  r  L,a+l,n+l 

w0  -  - I  W0  +  (  C1W1  +  C2w2  +  C3w3  +  C4w4 

C0T 

♦  S-+1,n+1}/  Cn  ♦  (  0-E“+1,n+1  +  (1  -  O  )  E*’n+1>/  C  (4.10) 

0  2  2  0 

For  <  1,  the  differencing  of  the  correction  ter»  becaae  a  weighted 


i 


average  of  Causs-Siedel  (current  level)  and  Jacobi  (previous  level) 
iterations.  Cenerally,  . 1  <  5  .5  was  used  in  the  numerical 

computations.  For  further  details  on  the  computational  procedure,  see 
Fant  (1987). 

Computational  Details  and  Discussion.  All  computational  runs  were 
made  using  31  radial  nodes  and  241  angular  nodes  for  the  entire  annulus 
(0  -  2;:).  The  mesh  spacing  was  slightly  compressed  by  a  factor  of  1.5 
near  the  inner  and  outer  walls. 

a 

For  all  calcualtions,  a  scaled  time  step  (At)  of  0.5  was  chosen. 

This  size  time  step  translates  to  small  actual  time  increments  since 
*  i/2 

t  =  (Ra/Pr)  t.  For  each  test  case  (inner  cylinder  temperature 
distribution),  C  values  of  increasing  magnitude  were  considered  until 
unsteady  behavior  resulted.  All  test  cases  were  started  by  obtaining  the 

A  A 

G  =  1.0  flow  field.  Then  for  each  slightly  higher  value  of  G,  the  flow 
conditions  from  the  previous  G  value  served  as  reference  data  for  the  new 

A 

set  of  results.  When  flow  instability  developed  at  some  value  of  G,  the 
number  of  iterations  for  convergence  and  the  computer  CPU-time  increased 
significantly.  Only  one  test  case,  Tj(*0  =  sin  y,  was  used  to  gain 
unsteady  flow  results  (see  Chapter  5).  All  other  test  cases  were  studied 
within  the  steady-state  regime. 

Generally,  steady-state  solutions  were  achieved  within  40  time 
steps,  although  solutions  were  taken  to  500  time  steps  to  ensure 
convergence.  The  amount  of  CPU-time  required  to  achieve  the  steady-state 
condition  was  approximately  38.4  seconds,  or  an  average  of  0.96  seconds 
per  time  step  on  the  AFIT  VAX  8650  computer  system.  Conversely,  about  4 
hours  CPU-time  was  required  to  run  200  time  steps  when  dealing  with  the 


unsteady  Multicellular  flow  regiae.  Such  tlae-consuiing  operations 
restricted  the  aaount  of  unsteady  flow  studies  that  could  be  conducted 
for  this  thesis. 
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V.  Results  and  Discussion 

In  this  chapter,  the  flow  fields  resulting  from  various  temperature 
distributions  imposed  on  the  inner  cylinder  wall  are  examined.  For  each 
variable  inner  wall  temperature  distribution,  Tj(v'),  a  set  of  isotherm 
and  streamline  contour  plots  are  presented  for  flow  visualization.  In 
the  streamline  contour  plots,  the  numbers  embedded  in  the  streamlines 
denote  stream  function  strength,  where  negative  values  denote 
counter-clockwise  fluid  rotation  and  positive  values  denote  clockwise 
rotation.  Positive  and  negative  isotherm  lines  represent  hot  and  cold 
regions,  respectively,  on  the  inner  cylinder.  Steady-state  analytical 
results  are  directly  compared  using  plots  of  vorticity  and  stream 
function  data  in  areas  of  interest  within  the  annular  flow  field.  In 
addition,  the  unsteady  flow  regime  shall  be  examined  for  one  of  the  test 
cases  described  below. 

Inner  Cylinder  Temperature  Distributions 

Figure  5.1  displays  the  ten  inner  cylinder  temperature  distributions 
considered  in  this  study.  Test  cases  (a  -  f)  were  chosen  to  impose 
regions  of  gradual  temperature  variation  around  the  annulus,  as  modeled 
by  simple  trignometric  functions.  More  localized  hot  and  cold  spots  were 
modeled  in  cases  (g  -  j).  The  presence  of  localized  heating  on  the  inner 
cylinder  wall  reflects  a  more  realistic  condition  than  the  isothermal 
boundary  condition  that  dominates  similar  studies  in  this  area.  Note 
that  the  circles  in  Figure  5.1  represent  the  inner  cylinder  with  the 
important  features  of  the  temperature  distribution  drawn  in  for  clarity. 
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Figure  5.1.  Inner  cylinder  temperature  distributions 
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The  functional  fora  of  the  annular  tenperature  variation  is  presented  to 
the  left  of  sketch.  More  detailed  isotherm  representations  of  these  test 
cases  are  given  in  the  next  section  and  in  Appendix  A. 

Humeri  cal  Results 

In  this  section,  several  test  cases  which  illustrate  the  variety  of 
flow  characteristics  encountered  in  this  study  are  discussed.  Results  of 
other  temperature  distributions  are  mentioned  only  briefly,  with  their 
stream  function  and  vorticity  plots  included  in  Appendix  A  for 
completeness.  In  all  but  a  few  cases,  numerical  results  were  obtained 
for  (scaled)  gap  numbers  of  3.0  and  4.40  to  examine  low  and  moderate 
buoyancy  effects  on  the  flow  circulation  within  the  annulus.  Three  inner 
wall  temperature  distributions  were  examined  up  to  the  unsteady  flow 
transition  point,  with  one  case  ran  at  a  gap  number  high  enough  to  induce 
unsteady  effects. 

Steady-State  Results.  Case  a,  Tj  «  sin  y,  was  found  to  yield 
steady-state  results  for  gap  numbers  up  to  4.95.  The  isotherms  for  this 
temperature  distribution  are  shown  in  Figure  5.2a.  For  gap  numbers  up  to 
4.93,  the  flow  field  exhibited  the  familiar  kidney-shaped  cell  pattern, 
shown  in  Figure  5.2b,  associated  with  an  isothermal  inner  cylinder 
(Tj  =  l).  As  Figure  5.2c  illustrates,  a  steady  multicellular  flow 
develops  in  the  vertical  portions  of  the  annulus  when  gap  number  is 
increased  to  4.95.  This  type  of  flow  instability  appears  to  be 
hydrodynamic  in  origin,  induced  by  fluid  buoyancy  shear  effects.  The 
instability  occurs  near  the  hottest/coldest  parts  of  the  inner  cylinder 
(y  *  90°,  270°).  This  makes  sense  physically  since  this  is  where  the 
largest  temperature  difference  between  the  inner  and  outer  cylinders 
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Figure  5.2J.  Angular  variation  of  Uie  center-line  stream  function 
for  T.(V')  =  sin  y/  at  G  =  4.95 


occurs  for  this  case  (see  Figure  5.2a),  thus  providing  the  maximum 
buoyancy  force  for  flow  convection  in  the  annulus.  Note  that  both  cells 
are  rotating  in  the  same  manner,  clockwise.  Unsteady  results  for  this 
case  using  C  =  4.99  are  discussed  in  section  5.2.3. 

Kidney-shaped  cellular  patterns  were  also  prevalent  in  cases 
c  and  d,  where  the  upper  half  of  the  inner  cylinder  was  cooled  or  heated 
while  maintaining  isothermal  conditions  (Tj  =  1)  on  the  lower  half. 

Here, the  cosine  function  wa3  employed  to  produce  the  greatest 
nondimensional  temperature  diference  between  the  inner  and  outer 
cylinders  at  y  =  180°,  where  Tj(rc)  =  0  for  the  cooled  case  (case  d)  and 
T j ( f7 )  =  2  for  the  heated  case  (case  c).  As  shown  in  Figures  ( A 1  a  -  c) 
the  effect  of  cooling  on  the  upper  half-cylinder  was  to  shift  the 
convective  circulation  toward  the  bottom  of  the  annulus.  The  opposite 
effect  occurs  when  the  upper  half-cylinder  is  heated  (Figures  A2  a  -  c). 
Increasing  G  from  3.0,  increased  the  stream  function  strength  but  did  not 
change  the  basic  cell  structure.  Another  type  of  flow  field  encountered 
was  the  formation  of  discrete  counter-rotating  cells  induced  by 
alternating  hot  and  cold  regions  distributed  evenly  around  the  inner 
annulus  wall.  The  temperature  distribution  of  Tj  =  -  cos(2-stO  shown  in 
Figure  5.3a,  produced  the  most  complex  of  these  multicellular  flow 
patterns.  Figures  5.3  b  and  c  show  the  six-cell  streamline  contours  for 
gap  numbers  of  3.0  and  4.40,  respectively.  Although  the  temperature 
peaks,  which  occur  every  90  degrees,  produce  identical  nondimensional 
inner  and  outer  wall  temperature  differences  of  1,  the  convective  cells 
at  ip  -  90°  and  270°  are  much  larger  than  the  other  four.  This  phenomenon 
is  due  to  the  influence  of  buoyancy  effects  being  strongest  at  these 
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annular  locations,  as  explicitly  shown  by  the  analytical  result: 

F1  or  sin  yi  (ly*?))  . 

Also,  vertical  symmetry  of  the  cell  formation  is  recognized. 

Other  test  cases  exhibiting  this  flow  behavior  were  the  inner  wall 
temperature  distributions  Tj  =  -  cos  y>  and  Tj  =  cos  y,  which  produced 
four  counter-rotating  cells  of  equal  size  and  strength.  No  cells  were 
predominantly  larger  due  to  the  equal  buoyancy  influence  at  each  cellular 
location  in  the  annulus.  Unlike  sin  yj ,  the  cos  y>  temperature 
distribution  did  not  induce  a  more  complex  multicellular  flow  field  when 

a 

the  unsteady  flow  transition  point  was  reached  at  C  =  5.75.  The  isotherm 
and  streamline  contour  plots  for  =  cos  y>  and  Tj  =  -  cos  yj  are  found  in 
Appendix  A  (Figures  A3  a  -  c  and  A4  a  -  c). 

The  test  cases  involving  more  localized  temperature  changes  on  an 
otherwise  isothermal  inner  cylinder  (cases  g  -  J),  model  hot  and  cold 
spots  near  the  top  of  the  annulus.  Case  i  models  a  triple  hot  spot 
distribution  symmetric  about  y>  =  180°,  where  a  temperature  peak  of  2.5  is 
reached.  Two  secondary  peaks  of  1.8  are  located  ±  22°  from  the  vertical 

symmetry  line  (see  Figure  5. A  a).  As  shown  in  Figures  5. A  b  and  c,  these 

hot  spots  induce  the  beginnings  of  a  secondary  cell  formation  near  the 
top  of  the  annulus.  This  flow  field  is  similar  to  the  isothermal  case, 
except  for  this  localized  increased  convective  activity. 

The  cold  spot  of  case  J  (Figure  5.5  a)  gives  rise  to  a  thermal 
instability  of  two  discrete  counter-rotating  cells  near  the  top  of  the 
annulus.  As  gap  number  is  increased  from  3.0  to  A. AO,  these  small  cells 

become  stronger  as  does  the  larger  kidney-shaped  cells  which  dominate  the 
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Figure  5.5 
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Figure  5 ,5e.  Angular  variation  of  the  center-line  stream  function  at 
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flow  field  (see  Figures  5.5  b  and  c).  Note  that  this  convective  behavior 
resembles  the  thermal  instability  of  air  ( Pr  =  .7)  in  an  isothermally 
heated  narrow  annulus. 

Cases  g  and  h,  presented  in  Figures  A5  a  -  c  and  A6  a  -  c  in 
Appendix  A,  model  single  and  dual  hot  spots,  respectively.  These 
temperature  distributions  produce  simply  an  elongation  of  the  upper 
portion  of  the  kidney-shaped  counter-rotating  cells.  Strong  secondary 
cell  formation  does  not  occur  as  in  case  i,  where  the  sharp  central  peak 
contributes  to  this  condition.  As  with  cases  i  and  j,  increasing  gap 
number  to  A. 40  only  increases  the  stream  function  strength  without 
altering  the  basic  cell  pattern. 

Analytical  Comparison .  The  analytical  perturbation  expressions 
derived  in  Chapter  III,  from  the  zero-Prandtl  number  boundary-layer 
equations  were  used  as  checks  to  the  steady-state  numerical  results. 

These  comparisons  support  both  the  flow  fields  found  by  the  numerical 
procedure  and  the  analytical  analysis  itself. 

The  two-term  perturbation  expressions  for  stream  function  and 
vorticity  were  used  to  generate  comparative  data  for  the  pretransi t ional 
numerical  results.  Stream  function  and  vorticity  were  compared  at  the 
angular  coordinate  where  maximum  stream  function  occurred  for  0  <  r  S  1. 
Comparison  of  stream  function  data  was  also  made  at  the  annulus 
centerline  (r  =  .5)  for  180°  <  v'  -  360°.  For  each  test  case,  the 
comparative  data  was  generated  for  the  largest  gap  number  used  in  the 
numerical  analysis  of  section  5.2.1. 

The  temperature  distribution,  T^CvO  =  sin  v,  demonstrates  the  limits 
of  the  analytical  formulation  in  describing  pretransitional  flows. 


Figure  5.2  d  shows  how  the  analytically  generated  streamlines  for 
G  =  4.95  resemble  the  numerical  result  of  G  =  4.93  (Figure  5.2  b),  while 
failing  to  capture  the  numerically-generated  multicellular  flow  structure 
of  Figure  5.2  c  (for  G  =  4.95).  Figures  5.2  e  -  g  explicitly  show  the 

A 

analytical  agreement  with  C  =  4.93  numerical  data,  while  Figures  5.2  h  -j 
highlight  the  descrepancies  in  the  G  =  4.95  results.  analytical  vs. 
numerical  comparison  plots  for  stream  function  and  vorticity  (Figures  h 
and  i)  indicate  significant  differences  in  the  solution  where  stream 
function  strength  is  maximum.  This  difference  is  also  seen  along  the 
centerline  (Figure  5.2  j)  where  descrepancies  exist  for  225°  <  y/  <  300°. 
This  is  the  multicellular  region  which  could  not  be  resolved  by  the 
two-term  perturbation  result. 

All  other  cases  exhibited  very  close  agreement  when  comparing 
analytical  and  numerical  solutions,  except  for  case  c;  the  cooled  upper 
half-cyl inder .  As  Figures  A1  d  -  g  show,  the  descrepancies  were  not  as 
dramatic  as  the  sin  w  case,  particularly  for  the  centerline  comparison. 
These  descrepancies  were  due  to  the  use  of  the  transitional  gap  number 
for  this  case  (G  =  4.71)  to  generate  comparison  data.  Although  the  cos 
test  case  was  evaluated  at  its  transitional  gap  number  (as  sin  v  was  for 
G  =  4.95),  no  discernable  differences  in  the  data  was  noted.  From  this, 
one  might  conclude  that  discrete,  counter-rotating  cells  tend  to 
stabilize  the  flow  field  compared  to  the  kidney-shaped  cell  pattern, 
since  the  comparison  data  for  the  other  seven  cases  were  within  the 
steady-state  or  pretransitional  zone,  no  descrepancies  were  found  (as 
expected).  The  comparison  plot3  are  included  with  each  test  case 
following  the  flow  field  contour  plots. 
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Unsteady  Results .  Only  the  sin  y/  inner  cylinder  temperature 
distribution  was  numerically  studied  in  the  unsteady  flow  regime.  As 
mentioned  in  section  5.2.1,  a  gap  number  of  4.99  induced  unsteady  flow 
behavior  in  the  annulus.  The  solution  was  carried  out  to  t  =  1000  to 
further  study  this  unsteady  development.  Flow  field  plots,  corresponding 

A 

to  t  =  200,  400,  600  and  800  are  presented  in  Figures  5.6  a  -  d.  As  time 
progressed,  the  weaker  cells  diminished  and  softened  in  structure.  Then, 
as  seen  in  Figure  5.6d,  some  re-strengthening  of  the  cell  pattern 
occurred  near  t  =  800.  Due  to  time  considerations  (on  the  average,  200 
time  steps  required  4  hours  of  CPU  time),  the  solution  was  not  carried 
out  further  to  determine  if  this  cyclic-like  flow  behavior  persisted. 
Although,  it  appears  that  some  type  of  cellular  structure  interaction  is 
occurring  in  the  vertical  portion  of  the  annulus. 
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VI .  Conclusions  and  Recommendations 

The  Influence  of  variable  inner  wall  teaperature  distributions  have 
direct  effects  on  the  formation  of  multicellular  convection  flows  in 
|  r-  narrow  horizontal  concentric  annuli  for  both  steady-state  and  unsteady 

flow  regimes.  For  most  test  cases  examined  in  this  study,  the  two-term 
analytical  perturbation  expressions  of  Chapter  III  provided  excellent 
j  agreement  with  pretransitional  numerical  results  for  steady-state  flows. 

These  solutions  begin  to  diverge,  however,  when  the  point  of  transition 
to  the  unsteady  flow  regime  is  reached,  as  shown  explicitly  by  the  sin  y/ 

|  '  test  case  of  chapter  V.  This  is  to  be  expected,  since  the  analytical 

formulation  is  derived  on  the  assumption  of  steady-state  conditions  and 
cannot  capture  any  time-dependent  multicellular  flow  field  behavior. 

|  |  These  solutions  are  therefore  best  suited  for  quick  flow  field 

visualization  of  steady-state  convection  in  the  annulus  subject  to 
variable  inner  wall  temperature  conditions. 

|  1  Abrupt  temperature  changes  on  the  upper  cylinder  surface  (hot  or 

cold  spots)  give  rise  to  small,  discrete  convective  cells,  especially  in 
the  horizontal  section  of  the  annulus,  which  resemble  thermal 
instabilities  characteristic  of  air  in  a  narrow  horizontal  enclosure. 

This  is  in  contrast  to  the  instability  formed  in  the  vertical  portions  of 
the  annulus  by  isothermal  or  gradually-varying  inner  wall  temperature 
conditions.  This  instability  is  hydrodynamic  in  origin,  induced  by 
buoyancy  forces. 

In  the  context  of  a  simplified  set  of  governing  equations,  one  can 
obtain  a  multitude  of  cellular  structures  in  the  annular  flow  field. 
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Such  complex  convective  flow  behavior  *ay  enhance  heat  transfer  of 
systems  eaploying  the  concentric  cylinder  geometry  for  cooling  or 
insulation  purposes.  Understanding  the  inner  wall  conditions  that  cause 
such  phenoaenon  is  central  to  this  analysis. 

In  order  to  gain  further  insight  into  the  formation  of  discrete 
aulticellular  steady-state  flow  fields,  more  studies  are  recommended  in 
the  pretransitional  regime  using  a  variety  of  hot  and  cold  spots  on  the 
lower  and  side  portions  of  the  inner  cylinder.  In  addition,  more 
numerical  work  should  be  performed  in  the  unsteady  and  transitional  flow 
regimes  for  a  variety  of  inner  cylinder  temperature  distributions  to 
study  the  nature  of  unsteady  multicellular  instabilities.  Experimental 
studies  of  low  Prandtl -number  fluid  convection  are  strongly  recommended 
for  the  narrow-gap  concentric  cylinder  geometry  with  variable  inner  wall 
temperature  conditions  to  verify  the  limited  numerical  work  performed  in 
this  area. 

Additional  numerical  and  analytical  work  is  also  recommended  to 
study  how  finite  Prandtl -number  fluids  effect  the  solutions  obtained  for 
the  test  cases  examined  in  this  thesis.  Such  research  would  offer 
insight  into  the  behavior  of  viscous  fluids  under  variable  inner  wall 
temperature  boundary  conditions  within  a  narrow  annular  enclosure. 

Finally,  a  parametric  study  of  the  numerical  grid  could  be  performed 
to  see  how  grid  density  affects  the  flow  field  solutions,  perhaps 
indicating  an  optimal  configuration  for  resolving  multicellular  flow 
behavior. 
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Appendix  A 


Supplementary  Flow  Field  Data 
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Figure 
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Figure 


numerical 


Figure  Ale.  Angular  variation  of  the  center-line  stream  function  at  ( 
for  Tjiy/)  =  l  +  cos  if/  imposed  on  the  upper  half-cylinder 


Figure  Alf.  Radial  variation  of  vorticity  at  88°  for  G  =4.71  for 
Tj {<(/)  =  1  ♦  cos  iff  imposed  on  the  upper  half-cylinder 
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Figure  A2e.  Angular  variation  of  the  center-line  stream  function  at  ( 
for  T j(y/)  =  1  -  cos  imposed  on  the  upper  half-cylinder 


nunerical 


[ill] 


Con  of  vorticity  at  120°  for  G  =  4.40  for 
os  */  imposed  on  the  upper  half-cylinder 


oro 


Figure  A2 
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numerical 


Figure  A3f.  Radial  variation  of  vortlcity  at  AS  *  for  G  =  5.75 
T.  (  r')  =  cos  ( V-') 
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Figure  A4e.  Angular  variation  of  the  center-line  stream  functi 
T.(^)  =  -  cos(V') 
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Figure  A5f. 
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Abstract 


</ 

The  purpose  of  this  study  was  to  examine  natural  convection  within  a 
narrow  horizontal  annulus  subject  to  variable  inner  cylinder  temperature 
distributions.  Both  numerical  and  analytical  approaches  were  taken  in 
determining  the  effects  of  variable  temperature  conditions  imposed  on  the 
inner  cylinder  in  triggering  or  suppressing  multicellular  flow 
instabilities. 

The  two-dimensional  Navier-Stokes  equations  are  simplified  into 
boundary-layer  equations  for  the  assumptions  of  large  Rayleigh  number, 
small  annular  gap,  and  small  Prandtl  number.  These  2-D  unsteady 
boundary- layer  equations  are  discretized  using  finite-differencing 
techniques.  Numerical  solutions  to  these  governing  equations  are  then 
obtained  by  using  a  stable  second-order,  fully-implicit ,  time-accurate, 
Causs-Siedel  iterative  procedure.  In  addition,  analytical  steady-state 
solutions  to  these  simplified  equations  are  obtained  using  perturbation 


For  most  inner  cylinder  temperature  distributions  considered  in  the 
steady  flow  regime,  these  analytical  results  yield  excellent  agreement 
with  numerical  results.  Although  both  schemes  predict  the  formation  of 
thermal-like  instabi 1 i ties  induced  by  localized  hot  and  cold  spots  in  the 
horizontal  annular  regions,  the  analytical  model  failed  to  predict  the 
steady-state  multicellular  hydrodynamic  instabilities  found  numerically 
for  the  sin  y  temperature  distribution  at  G  =  4.95  in  the  vertical 
portions  of  the  annulus.  The  analytical  model  also  fails  to  capture 
unsteady  multicellullar  flow  behavior  found  numerically  for  the 
sinusoidal  temperature  distribution  at  G  =  4.99. 


I 


